Fluctuations and topological transitions of quantum Hall stripes: 
nematics as anisotropic hexatics 
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We study fluctuations and topological melting transitions of quantum Hall stripes near half-filling 
of intermediate Landau levels. Taking the stripe state to be an anisotropic Wigner crystal (AWC) 
allows us to identify the quantum Hall nematic state conjectured in previous studies of the 2D 
electron gas as an anisotropic hexatic. The transition temperature from the AWC to the quantum 
Hall nematic state is explicitly calculated, and a tentative phase diagram for the 2D electron gas 
near half-filling is suggested. 
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Introduction - Following theoretical predictions by 
Koulakov et al. [H that the ground state of the 2D elec- 
tron gas near half filling of intermediate Landau levels 
(LLs), with index > 2, is a striped state, and the 
subsequent experimental observation by Lilly et al.j^ of 
strongly anisotropic dc resistivities in the above men- 
tioned range of fillings, it has been suggested [3, Q that 
the striped ground state of a two-dimensional electron gas 
(2DEG) at low temperature may be viewed as a "quan- 
tum Hall smectic" (QHS), consisting of a weakly cou- 
pled stack of one-dimensional Luttinger liquids. This is a 
state that would only be stable at zero temperature, and 
which, by analogy with conventional liquid crystals 
would give way through the proliferation of dislocations 
(see panels (a) and (b) of Fig. [T|) to a "nematic" state at 
nonzero temperatures Q , in which translational symme- 
try is restored but rotational symmetry is still broken. 
This electronic "nematic" would then undergo a discli- 
nation unbinding transition into a fully isotropic fluid as 
temperature is raised above a critical temperature which 
has been estimated 0, @] following standard Kosterlitz- 
Thouless (KT) arguments Q. 

In this paper, we want to examine an alternative pic- 
ture, in which the ground state of the 2DEG near half 
filling of intermediate LLs is taken to be an anisotropic 
Wigner crystal (AWC), as suggested by Hartree-Fock 
(HF) 1, [lo| and renormalization group (RG) [TH cal- 
culations. In this case, we find that dislocations melt the 
AWC at a nonzero temperature (that we shall explicitly 
evaluate below) into a "nematic" state with quasi-long- 
range orientational correlations, which we argue is noth- 
ing more than an anisotropic hexatic. Our results for the 
melting temperature of the AWC are consistent with ex- 
periments and with the idea of quantum Hall "nematics" 
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FIG. 1: Schematic view of (a) a smectic; (b) a smectic with 
an edge dislocation; and (c) an anisotropic Wigner crystal 
where the guiding centers of the electrons occupy the sites of 
an anisotropic triangular lattice. 



conjectured in Refs. 

Fluctuations of quantum Hall Wigner crystals ~ In 
what follows, we shall be interested in the elastic fluc- 
tuations of quantum Hall Wigner crystals. To fix ideas, 
we shall focus on the AWC which was found to mini- 
mize the cohesive energy of the 2DEG near half filling 
if Ref. 10, and which is described by the lattice vectors 
-n2a2, where ai = 2Q:y and Sl2 = ay-f /3x, 
aVl — e/2 and /? = -\/3a/2-\/l — e (ni and 
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^nin2 — ^l^r 

and where a = 



77.2 being integers). In these expressions, e is a positive 
parameter such that < e < 1 which quantifies the de- 
gree of anisotropy of the lattice at a given partial filling 
factor V* , and which was determined through minimiza- 
tion of the cohesive energy of the system in Ref. [l^; and 
a = £(47r/-\/3r^*) is the average spacing of a hexagonal 
lattice with e = at the same value of v* . The elastic 
properties of such an anisotropic crystal can be described 
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by an elastic Hamiltonian of the form (a,/? = x, y): 

J drdr' Cai3js{r - r')uaf3ir)ujsir'), (1) 

where Ma^(r) = ^{daUfj + dpUa) is the hnear strain ten- 
sor, (u(r) being the displacement field). For the par- 
ticular case of a two-dimensional AWC, there are three 
compression moduli, cn = Cim, C22 = C'2222 and 
C12 = C1122, and a single shear modulus Cee = (^1212 ■ 

The elastic fluctuations of the above AWC, taking into 
account the Lorentz-force dynamics imposed by the ex- 
ternal magnetic field, can be described by the Gaussian 
action: 

where the dynamical matrix _DQ,^(q, Wn) is given by: 

-Da/3(q, ^n) = <i>a/3(q) + Pm^^n^af3 + SapPrn^^c^Jn- (3) 

In the above expression, w„ = 2'imkBT /h is a Matsubara 
frequency (fcs being Boltzmann's constant and T being 
temperature), pm is a mass density, is the cyclotron 
frequency, Sap is the two-dimensional version of the an- 
tisymmetric Levi-Civita tensor, and we have introduced 
the elastic matrix $ct/3(q), which has the following matrix 
elements: 



J2 , „ „2 



(4a) 

$aa(q) = C22g^ +C66(?^, (4b) 
$:rj;(q) = $j^i(q) = (ci2 + cee)qxqy (4c) 



From Eq. we can easily derive the follow- 

ing expression for the two-point correlation function 

{ua{(l,ujn)up{q' ,u;i)) = {2nfS„.-iS{q + q')hGap{q,i^n), 

(5) 

where Sn.i is the Kronecker symbol, and where the prop- 
agator Gaf3 has the following elements (we sum over re- 
peated indices): 



Gapiq, UJn) 



with 



eQ7£/3i5'&7(5(q) + Pm^n^aP — SaPPm^c^-n 



Dei{D) 



(6) 



Det(L') = p^w^ + p^Lol [pmujl + $x4q) + *aa(q)] 

+ <i>Uci)%y{ci)-^ly{<D. (7) 

In real space, the mean squared displacement 
(MQ,(r, r)M/3(r, t)) can be written in the form 
(M„(r,r)M^(r,T)) = keT J^Gap{q), with the ef- 
fective propagator Gapiq) = Z^JT^-oo ^a/3(q' '^n)- 
Performing the Matsubara sum, we find that the static 



0.0050 
0.0025 
0.0000 
-0.0025 
-0.0050 
-0.0075 
-0.0100 



♦ ♦ 



♦ ♦ ♦ ♦ ♦ i' 



■ ci.i 



0.42 



0.43 



0.44 



0.45 



0.46 



0.47 



0.48 



FIG. 2; (Color online) Finite parts Cij, i, j = 1, 2, of the 
compression moduli for the anisotropic Wigner crystal as a 
function of the partial filling factor v* in Landau level N = 2. 



fluctuations of the quantum system at finite temper- 
atures can be described using the partition function 
Zci — /[(iu(r)] e"^/*"'^-^, with the effective classical 
Hamiltonian: 



H 



1 



i„(q)[G ^q)] upi-q). 



(8) 



Knowledge of the effective propagator G(q) allows us to 
study the effect of Lorentz-force dynamics on the elastic 
properties (and hence on possible topological transitions) 
of the system. Calculating the inverse effective propaga- 
tor G^^(q) and expanding the resulting expression near 
g = 0, we find that, up to terms of order 0{q^), the form 
of the elastic propagator is identical to its zero field ex- 
pression. This has the important consequence that the 
long wavelength elastic properties of the AWC will be 
qualitatively the same as in the absence of magnetic field. 
We therefore expect the topological melting of the AWC 
to proceed in a standard (two-stage) way [16], as we now 
are going to describe. 

Topological melting of anisotropic Wigner crystals - 
A major difference between isotropic and anisotropic 
Wigner crystals in two dimensions is that, while in the 
former all six elementary dislocations differing by the 
orientation of their Burger's vectors are equivalent, in 
the latter two elementary equivalent dislocations (labeled 
type I) have their Burger's vectors along a reflection sym- 
metry axis (i.e. along ±ai in Fig. 1), while four dislo- 
cations, equivalent to each other but inequivalent to the 
first type, lie at angles of ±6*0 from the reflexion axis {9q 
here is the angle between ai and a2, see Fig. [IJ. At any 
nonzero temperature, the solid phase has a finite density 
of tightly bound dislocation pairs. As the temperature is 
raised past a critical temperature Tci, the pairs unbind 
and destroy the crystalline order. Since the two types 
I and II of dislocations are unequivalent, the defect me- 
diated melting (DMM) process will be governed by the 
type which has a lower nucleation energy. We therefore 
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FIG. 3: (Color online) Critical temperature Tci for the dislo- 
cation mediated melting of the anisotropic Wigner crystal into 
a nematic-like state near half filling of Landau level N = 2. 



shall need to determine the elastic constants of the AWC 
in order to find the energies of the two types of disloca- 
tions, so as to determine which dislocation type unbinds 
first. 

For the AWC which is the object of study in this paper, 
the compression moduli Cij{q), i.j = 1, 2, can be written 
in the form Cy(q) = c{q) + Cij, — 1,2), where we 
separated out the leading (plasmonic) contribution [15 1 



-/2TTKt j /q£, e being the electronic charge, 



the magnetic length, and k the dielectric constant of the 
host medium. For a one-dimensional compression of the 
form Ui = uqxSc, with the number of electrons in 
Landau level N kept fixed, if we denote by ;/* = ly* + 
uq) the partial filling factor of the compressed crystal, it 
can be shown [3] that the constant part cii(z^*) is given 
by: 



5ii {i^n 



{ev* 



dv* 



dv* 



(9) 



In the above expression, Q{v) is the HF cohesive energy 
per electron (in units of / nt), and is given by: 



(10) 



Q 



where Q is a reciprocal lattice vector, and p (Q) is the 
guiding-center density operator, which is determined self- 
consistently using the approach of Ref. 0, and is related 
to the real density operator n(Q) through the equation 
(here is the Landau level degeneracy and (x) is a 
generalized Laguerre polynomial): 



i(Q) = N^e-Q'''/^L% 



p{Q). (11) 



On the other hand, the Hartree and Fock interactions are 



given byQ (Jq (x) is the Bessel function of order zero): 



HN{q) 



le-9'^V2 



/o2^2 



(12a) 



Xn (q) = \/2 y dx e-^' [L% {x^)f Jq (^V2xq£ 

(12b) 

The finite contributions C22 and C12 to the compression 
moduli C22 and C12 can be obtained in a similar way by 
considering ID and 2D uniform compressions of the form 
u = uoyy and u = uo{xx + yy), respectively. The results 
of these procedures, the details of which will be published 
elsewhere IJ], are shown in Fig. ^ where we plot the 
constant parts ctj {i, j — 1, 2) of the compression moduli 
of the stripe crystal near half filling of LL N = 2. 

Let us now introduce the compliance tensor Sijki such 
that SijkiCkimn = ^{SimSjn + ^inSjm), with the four 
independent elements Sii(q) = ^nn, S22(q) = 52222, 
si2(q) = S2i(q) = 5'ii22, and sge = 51212. In the long 
wavelength limit, these take the values: 
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Sll = S22 



-S12 



cii + C22 - 2ci 



S66 



4cfi 



(13) 



In terms of these compliances, we find that the leading 
contribution to the energy of a dislocation of type a {a = 
I or //) is given by 



1 



K^HL/a), 



with hj: 




Kji -Ki 



where we defined 



K = — [2S22(\/S11S22 
ZTT 



S12 



2S66)] 



1/2 



(14) 

(15a) 
(15b) 

(16) 



For the problem at hand, sn = S22 in the long wave- 
length limit, and hence we see that the ratio Kjj/Kj is 
always larger than unity for < e < 1. We thus see that 
dislocations of type I are energetically less costly than 
type II dislocations, and will unbind at the melting tem- 
perature Tci such that KB\ Kjo? /ksTd = 4 (the value 4 
being universal), from which we obtain the melting tem- 
perature 



(17) 



The resulting melting line of the AWC is plotted in Fig. 
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FIG. 4: Tentative phase diagram showing the various relevant 
regimes for the melting of quantum Hall stripes as a function 
of temperature T and length scale L. ^disc is the average 
spacing between disclinations. Below this length scale above 
Tc2 the system retains the properties of an anisotropic hexatic. 
Note that Tci here can be extremely small, and may be driven 
to zero by quantum fluctuations. 

For temperatures higher than Td , the presence of type 
I dislocations screens the logarithmic interaction between 
type II dislocations, such that both dislocation types are 
free at long length scales [l^. We expect, however, that 
the average separation of type II dislocations will be 
much larger than the average separation ^/ between type 
I dislocations, with [l^: 

6 ~ exp(i-i/2), 6/-^r\ (18) 

where t on {T — Td) and p is a nonuniversal number 
between and 2. Given that type I dislocations de- 
stroy translational order mainly along the direction of 
the stripes, we see that we can distinguish between three 
different regimes (see Fig. At length scales shorter 
than the system retains the properties of a solid. At 
intermediate scales, ^/ < L < ^/j, the system is smectic- 
like, and consists of a regular stack of ID channels of 
electron guiding centers, with short ranged translational 
order along the channels, and quasi-long ranged order in 
the transverse direction. Finally, at length scales longer 
than ^/7, the system is nematic-like: translational order 
is destroyed in all directions, but the system preserves 
quasi long range orientational order. The latter is de- 
scribed by the bond-angle field 9{r), which is defined as 
the orientation relative to some fixed reference axis of the 
bond between two neighboring electron guiding centers. 
Standard analysis shows that the fluctuations of 9{r) are 
described by an effective Hamiltonian of the form: 

= i y rfr [K.,{d^0f + Ky{dy9f] , (19) 

and decay only algebraically with distance. Since on 
short length scales, the AWC is only weakly disturbed 
and each electron is still surrounded on average by six 



neighbors, he resulting quantum Hall "nematic" state 
may be more accurately characterized as an anisotropic 
hexatic. 

As the temperature is further raised, a disclination- 
unbinding transition melts this nematic-like state into 
an isotropic metallic state, in much the same way as de- 
scribed in Ref. |16|, with actual values of the nematic to 
isotropic melting temperature Tc2 of the order of those 
estimated in Ref. d (in this last reference, » 200 mK 
near half-filling of LL iV = 2). Since the temperature 
('^ 25mK) at which the experiments of Ref. Q were per- 
formed lies between Td and Tc2, we see that our HF 
calculation is consistent with the conjecture [1, |3| ac- 
cording to which the state probed by these experiments 
is a nematic state. 

Conclusion - To summarize, in this paper we have 
examined the fluctuations and topological transitions 
of quantum Hall stripes near half-filling of intermedi- 
ate Landau levels. Taking the stripe state to be an 
anisotropic Wigner crystal, as suggested by Hartree-Fock 
and renormalization group calculations, we find that the 
quantum Hall nematic conjectured in Refs. 0,0] emerges 
in a natural way in the topological melting process, and is 
identified as an anisotropic hexatic. Our calculations are 
consistent with the idea of quantum Hall nematics, which 
we predict to be realized over a significant region of the 
phase diagram near half filling of intermediate Landau 
levels, and give quantitative support to the qualitative 
interpretations [1, |3| of transport measurements f^] in 
terms of putative nematic states. 
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